A minisatellite-based MLVA for deciphering the global epidemiology of the bacterial cassava pathogen Xanthomonas phaseoli pv. manihotis

Cassava Bacterial Blight (CBB) is a destructive disease widely distributed in the different areas where this crop is grown. Populations studies have been performed at local and national scales revealing a geographical genetic structure with temporal variations. A global epidemiology analysis of its causal agent Xanthomonas phaseoli pv. manihotis (Xpm) is needed to better understand the expansion of the disease for improving the monitoring of CBB. We targeted new tandem repeat (TR) loci with large repeat units, i.e. minisatellites, that we multiplexed in a scheme of Multi-Locus Variable number of TR Analysis (MLVA-8). This genotyping scheme separated 31 multilocus haplotypes in three clusters of single-locus variants and a singleton within a worldwide collection of 93 Xpm strains isolated over a period of fifty years. The major MLVA-8 cluster 1 grouped strains originating from all countries, except the unique Chinese strain. On the contrary, all the Xpm strains genotyped using the previously developed MLVA-14 microsatellite scheme were separated as unique haplotypes. We further propose an MLVA-12 scheme which takes advantage of combining TR loci with different mutation rates: the eight minisatellites and four faster evolving microsatellite markers, for global epidemiological surveillance. This MLVA-12 scheme identified 78 haplotypes and separated most of the strains in groups of double-locus variants (DLV) supporting some phylogenetic relationships. DLV groups were subdivided into closely related clusters of strains most often sharing the same geographical origin and isolated over a short period, supporting epidemiological relationships. The main MLVA-12 DLV group#1 was composed by strains from South America and all the African strains. The MLVA-12 scheme combining both minisatellite and microsatellite loci with different discriminatory power is expected to increase the accuracy of the phylogenetic signal and to minimize the homoplasy effects. Further investigation of the global epidemiology of Xpm will be helpful for a better control of CBB worldwide.

Introduction Molecular genotyping of bacterial pathogens allows a comprehensive view of their epidemiology and evolution. The ability to differentiate genotypes at different spatial and temporal scales can be informative of the dispersion and micro or macroevolution of bacterial pathogens [1][2][3]. At a local scale, this approach can shed light on the microevolutionary processes within populations, and track back the spread of bacterial isolates and the associated pathways. At a long-term or global scale, it can for instance unravel the invasion routes and characterize some bacterial lineages associated to pathological or adaptative traits for epidemiological surveillance. Genotyping methods based on several molecular markers have been developed to analyze the genetic diversity of pathogenic bacteria. The type of marker should be selected according to the research question to be answered because each marker has different mutation rates i.e. different speed of evolution and stability [2][3][4][5][6]. Markers with a desired molecular clock can be selected according to the evolutionary scale under study.
Global epidemiology studies or international surveillance of pathogens require tools that can produce unambiguous data and are easily transferrable among different laboratories. Multilocus sequence typing (MLST) is a nucleotide sequence-based approach using markers that are conserved, stable, with relatively low evolution rates. It enables international comparison of isolates and global epidemiology [7]. However, some bacterial pathogens, the so-called genetically monomorphic bacteria [8], do not exhibit enough variability in their housekeeping genes making this approach insufficiently resolutive [8]. Various subtyping methods emerged with the advent of whole genome sequencing. In situations where large numbers of high-quality draft genomes can be obtained, single nucleotide polymorphism (SNP) typing is a powerful approach for epidemiological surveillance and phylogenetic analyses [8][9][10][11]. Technically simple and still robust systems are also necessary for routine epidemiological investigations in situations where an easy access to large sets of compete genomes is impossible to achieve. In such situations, Tandem Repeats (TRs) and Clustered Regularly Interspaced Short Palindromic Repeats (CRISPRs) have been mainly applied as typing tools with regard to their reliability, discriminatory power, transferability, ease of application and cost [12,13].
TRs are short motifs of DNA fragments repeated in tandem that can vary in number due to addition or deletion of motifs. Multiple loci variable number tandem repeat analysis (MLVA) targets multiple TR loci which allows a high discriminatory power if a sufficiently large number of TRs is targeted. Two classes have been proposed to separate these informative molecular markers: microsatellites with a repeat unit size < 10 bp and minisatellites with repeat units � 10 bp [14,15]. The length of TR arrays is typically determined by capillary electrophoresis from multiplex involving fluorochrome-labeled primers but it can be scored unambiguously and from agarose gels especially for large repeat motifs making them applicable at a lower throughput in laboratories where a capillary electrophoresis genotyper is not available. Several microsatellite schemes have been developed for epidemiological analyses of plant pathogenic xanthomonads [16]. Minisatellite-based typing schemes targeting large repeat motifs were also proposed as a portable approach for the global analysis of the molecular epidemiology of bacterial strains pathogenic to human or animals [15,17,18] or to plant crops [19].
The wide range of mutation rates offered by the different TR loci classes make them useful for examining genetic patterns at several evolutionary scales. While microsatellites are most suited for local epidemiology studies, minisatellites have been mostly used for global epidemiology, i.e. analysis of non-epidemiologically related isolates [17,19].
Mutational processes observed in VNTRs are mainly driven by DNA replication slippage (and to a lesser extent recombination) and the length or the number of the tandem repeat and the sequence of the repeat unit are factors that drive the mutation rate of TRs [14,20,21].
Cassava Bacteria Blight (CBB) is a very important disease caused by Xanthomonas phaseoli pv. manihotis (Xpm) [22,23]. This disease was first reported in Brazil in 1912 [24], and later noticed in Colombia in 1970 [25], Venezuela in 1971 and African and Asian countries since 1972 [26,27]. CBB might have been present in these countries for many years, but few reports described a bacterial leaf spot on cassava before the seventies. More recently, the disease was reported in South Pacific [28]. Population studies of the pathogen have been conducted in several Latin American countries, including Brazil, Venezuela [29,30] and a continuous monitoring of the populations has been performed in Colombia [31][32][33][34][35][36]. Xpm is transmitted at small scale by rain and wind within production areas. At larger scales between production areas and between different countries, contaminated plant material, i.e. stem cuttings, are responsible for bacterial dispersal and the carry-over to the following growing season [24].
Studies of the genetic diversity in Xpm populations have been mainly performed at a local level using Restriction Fragment Length Polymorphisms (RFLP) or Amplified Fragment Length Polymorphisms (AFLP), as well as whole genome sequencing [37]. In Colombia Xpm populations were genetically diverse and geographically structured [34,38]. Genetic variation was observed within fields with high levels of genotypic diversity and haplotype frequencies significantly differing over the years [33,38]. More recent studies in different regions of Colombia confirmed the high genetic diversity and the local and geographic structure of Xpm populations in these regions [35,36]. A geographic structure of Xpm populations was also described in Venezuela at a regional scale [29]. Similarly, a high population diversity has been observed in African countries [39].
A few large-scale analyses of the genetic diversity of Xpm populations have been performed. Analysis of a worldwide collection using RFLP led to the identification of clusters grouping strains that originated from different countries [26]. This study suggested that African strains belong to a single subclade of South American strains, whereas Asian strains belong to different groups. A recent study based on single nucleotide polymorphism analyses supports gene flow among Xpm strains representative of the worldwide diversity especially between Brazil and Colombia and between the South-American and the African continents [37]. A common ancestor among Brazil, Colombia and African strains was suggested as previously suspected [26].
Herein, we built on the versatility of TR markers in order to develop improved MLVA schemes with desirable characteristics for outbreak investigation or epidemiological surveillance of Xpm. An MLVA-scheme targeting microsatellites was recently proposed and evaluated on local populations [40]. Our aims in this study were (i) to develop new TR markers with larger repeat size, (ii) to evaluate the discriminatory power of these minisatellites on a worldwide collection of Xpm strains in comparison with the MLVA scheme based on microsatellites and (iii) to propose an MLVA scheme adapted for regional or international surveillance of Xpm.

Xpm strain collection and DNA extraction
A collection of 93 Xpm strains isolated over several years from different cassava producing countries from South America (Colombia, Venezuela, Brazil, Argentina), Africa (Burkina Faso, Mali, Togo), Asia (Vietnam, and China) and New Zealand, (S1 Table) was used for genotyping. The genomic DNA was extracted from 24 hour-old cultures on an LPGA culture medium (Yeast extract 5 g L -1 , Peptone 5 g L -1 , glucose 5 g L -1 , agar 15 g L -1 , pH 7.2). DNA extraction was performed using the Wizard 1 Genomic kit (Promega, Madison, WI, USA). The amount and purity of DNA was checked using NanoDrop technology (Thermoscientific, Illkirch, France).

Genotyping using a new MLVA scheme targeting minisatellites
Twelve Xpm genomes, previously sequenced with Illumina [37], were selected, based on degree of fragmentation (least fragmented were included) and number of contigs (lower number of contigs were included). Eight strains from Brazil, two from Colombia and two from Uganda [40], were screened for VNTR minisatellite loci using the Xanthomonas utility website (http:// bioinfo-web.mpl.ird.fr/xantho/utils/), setting the following parameters: 30 and 1000 bp for the total length, the unit length of the repeat in a range of 10 to 100 bp, similarity of repeats within the arrays between 80 and 100%.
The primers were designed over flanking regions up-and downstream of the repeat region as previously described [40]. Primers were tested by individual PCR and multiplex PCR using DNA from a sub-collection of Xpm strains: CIO1, CIO151, COL303, CFBP1851 (Colombia), VEN130 (Venezuela), SCI (Ivory Coast), D2-3 (Burkina Faso), Mali45, Mali30 (Mali). These strains were only used to test the quality and design of the primers in agarose gels.
Multiplex PCR was performed to amplify three loci in each multiplex using the QIAGEN 1 Multiplex PCR kit (Qiagen, Courtaboeuf, France) according to the manufacturers´instructions with 40 to 50 ng of target in a total volume of 10 μl. The multiplex PCR conditions were as follows: pre-denaturation at 95˚C for 5 minutes, 30 cycles of denaturation at 95˚C for 30 seconds, annealing with a specific temperature per pool for 40 seconds (Table 1), extension at 72˚C for 2 minutes; and a final extension at 72˚C for 5 minutes. The amplicons were visualized by agarose gel electrophoresis (2%).
After the standardization of the multiplex PCR, the reverse primer in each primer pair was labelled using different fluorescent dyes to separate the PCR products by capillary electrophoresis (Table 1). PCR products were diluted 20 times (estimated from preliminary tests), and 1 μl of the dilution was mixed with HiDi formamide (9 μl) (Applied Biosystems, CA, USA), and GeneScan TM 600 LIZ 1 dye size standard v2.0 (0.5 μl) (Applied Biosystems, CA, USA). The detection was performed in an ABI 3500 Genetic Analyzer, and amplicon sizes were determined using GeneMapper1 (Applied Biosystems) software version 4.1. The strain CIO151 [41] was used as a reference strain in each test.

Xpm genotyping using MLVA-14 targeting microsatellites
The strain collection was also genotyped using a MLVA scheme targeting 14 microsatellites, referred to as MLVA-14. This scheme is derived from MLVA-15 [40], but a single microsatellite, Xpm 2-20 (VNTR Xpm 1-20 nomenclature from [40]), was taken out because it did not amplify on a few strains from Venezuela or Africa (unpublished data). One primer from each pair in the PCR mix was labelled with one of the four different fluorescent dyes (6-FAM, NED, PET, and VIC, Applied Biosystems) (S2 Table) and detected by capillary electrophoresis in an ABI 3500 Genetic Analyzer as indicated above.

Data scoring
Fragment sizes for each VNTR locus were estimated using GENEMAPER v 4.0 (Applied Biosystems). When necessary, a truncated number of repeats was rounded up to the nearest integer [12]. The discriminatory power of each VNTR locus was calculated based on the Hunter and Gaston discriminatory index (HGDI) [42]. Numbers of alleles and Nei's unbiased estimates of diversity were estimated using ARLEQUIN software [43]. Haplotype networks were represented by minimum spanning trees using the algorithm combining global optimal eBURST (goeBURST) and Euclidian distances in the software PHYLOViZ 2.0 [44]. The minimum spanning trees display the relationships between clonal complexes defined as groups of single locus variants (SLVs). The population structure of Xpm was assessed by discriminant analysis of principal components (DAPC) of MLVA-12 data using the adegenet V.2.1.1. R package [45,46]. DAPC is free of any assumption linked to a population genetic model (e.g., Hardy-Weinberg equilibrium or linkage equilibrium), and, thus, is suited for analysis of datasets produced from predominantly clonal bacteria.

Selection of eight new minisatellite loci
After the screening of 12 Xpm genomes, [49] VNTR loci met the criteria of the search. Of these candidates, only ten loci were selected because their flanking regions were almost fully conserved, and their sequences were identical for each locus in all the strains analyzed. Five TR loci did not show any polymorphism for these 12 genomes. For the polymorphic loci the HGDI was relatively high varying between 0.167 and 0.849 with an allele number ranging from two to six, respectively ( Table 2). The total length of most of the loci including the primer regions ranged between 197 to 490 bp. The unit length of the tandem repeat ranged between 10 and 26 bp. Individual bands were observed for each primer and in PCR multiplex (Fig 1). However, for loci Xpm2-29 and Xpm2-2, the size of amplicons exceeded the range size covered by capillary electrophoresis. The sequencing analysis of these two loci showed large insertions so they were removed from the typing scheme. The final minisatellite scheme was formed by eight VNTR loci and called hereafter MLVA-8 (Table 1). This scheme had three perfect loci (Xpm 2-20, Xpm 2-23, and Xpm 2-33) and five imperfect loci ( Table 2).

MLVA-8 genotyping of a world collection of Xpm strains
Thirty-one multilocus haplotypes were obtained from the worldwide collection of 93 Xpm strains using the MLVA-8 scheme. All the TR loci were polymorphic except Xpm2-5, which we hypothesize was fixated by natural selection. Polymorphic loci produced between two (loci Xpm2-33 and 2-35) and eight alleles (Xpm2-18) with Nei's gene diversity varying from 0.108 to 0.831 (S3 Table). TR markers along the evolutionary path of the minimum spanning tree were analyzed by assessing the number of repeats involved in the polymorphism of recently diverging haplotypes, i.e. SLVs. Multiple-repeat variations occurred only for loci Xpm 2-3 and  c Nomenclature from Trujillo et al. [35] and/or from Arrieta-Ortiz et al. [41] and/or Rache et al. [40]. https://doi.org/10.1371/journal.pone.0285491.t002

PLOS ONE
Xpm 2-20 with 50% of double-or triple-repeat variations and 33% of double repeat variations, respectively. Changes consisting of a single-repeat variation represented 85.2% of the SLVs. Our data suggest that these TR minisatellites mainly evolve following a generalized stepwise mutation model (or two-phase model) [48]. Eighteen haplotypes out of 31 originated from a unique country. Thirteen haplotypes were shared by strains originating from different countries, six of which were from South American countries and seven from different continents. All the African strains shared haplotypes with South American strains except MLVA-8 haplotype #5 composed only of Burkinabe strains. The minimum spanning tree separated the haplotypes in three clonal complexes (CC, i.e. clusters of single locus variants), and one singleton, i.e. haplotypes with no SLVs (Fig 2). The main MLVA-8 CC1 grouped 24 haplotypes including all the African strains and strains from the other countries included in the analysis, except a single strain from China. They were isolated from 1966 to 2016. Three MLVA-8 haplotypes were recovered over a period of more than 30 years (# 1, 7 and 16) (S1 Table and Fig 2). All haplotypes originating from Mali or from Burkina-Faso were SLVs (MLVA haplotypes #2 and 4 differ only by Xpm 2-18). The MLVA-8 CC2 (n = 4 haplotypes), a triple locus variant of CC1, was composed by the Chinese strain and strains from Brazil and New Zealand that were all isolated between 1974 and 1996. The third CC (2 haplotypes) and the singleton, haplotype #21, grouped strains from Colombia and New Zealand isolated between 1966 and 1974. Strains from China and Vietnam did not share any haplotypes with other strains.

MLVA-14 genotyping of a world collection of Xpm strains
We compared the discriminatory power of MLVA-8 with the previously reported MLVA-14 microsatellite scheme. MLVA-14 scheme was able to distinguish 89 haplotypes out of the 93 Xpm strains. All loci were highly polymorphic with a mean allelic richness of A = 13.5 leading to high discriminatory index (HGDI = 0.999) and Nei's gene diversity index (H E = 0.863), which were much higher than those estimated from the MLVA-8 dataset (Tables 3 and S4). Strains sharing the same haplotype originated from the same country and, when the information was available, they were isolated the same year in the same locality (S1 Table).
The minimum spanning tree describing the relationships between the 89 haplotypes produced ten small clonal complexes grouping only two or three SLVs (S1 Fig). All these CCs grouped epidemiologically related strains from the same region isolated the same year except

PLOS ONE
for the two MLVA-14 haplotypes #3 and 5 (CC7) from two different Colombian regions and isolated during two consecutive years (S1 Table). Other closely related haplotypes (i.e. DLVs) also originated from the same country with the exception of two strains from New Zealand isolated in 1966 (haplotypes #44 and 69) which are DLVs of strains from Venezuela and Colombia isolated in the beginning of the seventies, respectively (S2 Fig and S1 Table). Almost all the other haplotypes differed from each other by at least four TR loci but mostly gathered by   (Table 3). Seventy-eight MLVA-12 haplotypes out of 93 strains were separated into 16 CCs (n = 60 strains) composed of two to seven haplotypes and 26 singletons (n = 33 strains) (S2 Fig). Nine CCs contained strains originating from the same country among which six consisted of strains isolated the same year in the same region probably corresponding to the same outbreak. Five CCs grouped strains from the same continent either South America (CCs 6, 11 and 16, each with strains from two different countries) or Africa (CCs 1 and 2, both with strains originated from Mali, Burkina Faso and Togo). Most of these CCs are tightly related as they differ by variations in two loci. We used DAPC to delineate seven clusters in our dataset (Fig 3). The assignation probability to clusters was > 0.99 for all strains in the dataset. Interestingly, the large polymorphism revealed for strains from South America split them in the seven clusters. Strains originating from Argentina (n = 4), Brazil (n = 12), Colombia (n = 28) and Venezuela (n = 21) were assigned to three (Argentina) or four (other countries) distinct clusters. Cluster #1 included strains from South America (Argentina, Brazil and Colombia) together with strains from Viet Nam. Cluster #2 mostly included strains from South America (Argentina, Brazil and Venezuela) together with a single strain from New Zealand. All African strains gathered in cluster #3 together with some strains from Colombia and Venezuela. Clusters #4 and 5 solely included strains from Colombia and Venezuela, respectively. Cluster #6 primarily included strains from South America (Brazil and Colombia) together with a few strains from China and New Zealand. Cluster #7 included strains from all countries in South America (Argentina, Brazil, Colombia and Venezuela). The structure revealed by DAPC was overall consistent with the produced MST (Figs 3 and S2).

Pathogenicity patterns are not associated to genetic groups
The AUDPC values were used to determine if cassava varieties are resistant (∑AUPDC �39), moderately resistant (39<∑AUDPC<44), moderately susceptible (44<∑AUDPC<49), or susceptible (∑AUDPC�49) to the eight Xpm strains evaluated (Table 4 and S3 Fig). One cultivar, cv. 60444, was susceptible to all strains tested and two varieties (NGA11 and CM523-7) were susceptible to nearly all strains tested. On the other hand, cultivars COL1505 and CM6438-14 were resistant or moderately resistant to at least three of the strains evaluated (Table 4).
Regarding the strain pathogenicity level, strains CIAT1202, UA556, UA2164, CIO151 and CIAT1241 caused disease in all inoculated varieties. On the opposite hand, strain CIAT1135 was poorly pathogenic. Based on this information, pathological patterns were classified as follows. The more virulent strains (CIAT1202, UA556, UA2164, CIO151 and CIAT1241) were classified in group I. Strains with intermediate virulence (UA1591 and CIAT1205) were designated in group II. The least virulent strain CIAT1135 was assigned to group III (Table 4). Data shown is the mean of the sum of AUDPC in arbitrary units at 28 dpi from five replicates per genotype. A genotype was considered resistant (R) if ∑AUPDC �39, moderately resistant (MR) if 39<∑AUDPC<44, moderately susceptible (MS) when 44<∑AUDPC<49 and susceptible (S) when ∑AUDPC�49, NE: Not evaluated.

Discussion
Epidemiology can be deciphered at different time and space scales by selecting the level of genomic resolution of the markers to monitor either the geographic spread of inoculum or the prevalence of certain epidemic or endemic clones for surveillance [2][3][4]6]. Cassava Bacterial Blight is a destructive disease widely distributed in the different tropical areas where cassava is grown. It has recently been reported in new areas [50]. A global molecular epidemiology analysis is needed to better understand the expansion of the disease for improving the monitoring of CBB. We recently developed a MLVA scheme targeting microsatellites with a high discriminatory power expected to be more adapted to local epidemiology and outbreak analysis [40]. Here, we evaluated new minisatellites for their usefulness in describing the relationships of non-epidemiologically related strains and large-scale epidemiology. These tandem repeat sequences, i.e. micro-and minisatellites, show motifs whose length variability is one of the main factors known to influence the mutation rate of VNTRs [14]. We propose a MLVA-12 scheme which takes advantages of combining TR loci with different mutation rates for global epidemiological surveillance.
Eight minisatellites (MLVA-8) with a unit length ranging from 10 to 26 bp (i.e. markers with an expected low mutation rate) separated a world collection of Xpm strains isolated over a period of fifty years  in three CCs and a singleton. The major CC1 grouped strains originating from all the represented countries, except the unique Chinese strain that was included in our collection. All the strains from Africa grouped in this CC while strains from South America, the continent where the disease probably originated, were distributed in the three CCs. The MLVA-8 dataset confirmed (i) a more extensive genetic diversity for the South American strains than for the African strains and (ii) the hypothesis of some genetic links between African strains and some South-American strains as previously shown using four RFLP probes and SNPs based-phylogeny [26,37].
The MLVA-8 scheme alone lacked discriminatory power and therefore did not reveal a fully meaningful epidemiological perspective at a global scale. The relatively low number of minisatellite loci identified and their low mutation rates both are factors favoring homoplasy [51,52]. In contrast, 89 out of 93 Xpm strains genotyped using the MLVA-14 scheme were separated as unique haplotypes clearly confirming the high mutation rate of this genotyping technique, consistent with a previous study and suitable for small scale epidemiology [40].

PLOS ONE
A MLVA-12 scheme combining the eight minisatellites and four microsatellites from the MLVA-14 and displaying the lowest evolution rates is proposed to increase discriminatory power and minimize homoplasy effects [51,52]. In consequence the MLVA-12 scheme, will increase the genetic space of possible patterns as compared to MLVA-8 and would be suitable for epidemiological surveillance. The variable mutation rates within the MLVA-12 scheme allowed a nested approach where the genetic relationships among strains is described within deeper phylogenetic clusters. Such approaches, when combining markers with different discriminatory power, progressively or simultaneously, have been found to increase the accuracy of the phylogeny and to minimize the homoplasy effects [4,52,53].
A DAPC analysis of MLVA-12 data separated the haplotypes in seven distinct clusters. Interestingly, this analysis placed strains from South America in all clusters. This further supports the wide genetic diversity already described in South America (i.e. center of origin of cassava and the continent where CBB was firstly reported [26,30,33,37]. In contrast, all African strains studied here grouped in a single cluster also containing strains from South America, suggesting a possible relatedness and the existence of a common ancestor. However, based on MLVA-14, the amount of genetic relatedness between African and the closest South American strains did not support recent epidemiological links. Cassava was introduced in Africa during the seventeenth and eighteenth centuries at a period where regular communications already existed between the two continents [54]. After its first description in Brazil in 1912, cassava bacterial blight was first reported in continental Africa in the early 1970's [24,54]. The pathogenic bacterium was likely introduced in this continent through transportation of infected plant material after multiple introduction events into both Eastern and Western African coasts. Three strains isolated in New-Zealand during outbreaks in 1960s or in 1980 were distributed in two different DAPC clusters, one of them sharing an MLVA-12 haplotype with a strain from Venezuela and the others were SLVs of Colombian and Brazilian strains. This result supports the occurrence of multiple introduction events of Xpm in this country.
No clear correlation between the MLVA profiles and virulence patterns occurred. The different haplotypes representative of the genetic variability of observed within our world Xpm collection from MLV8-dataset produced different pathological patterns though all the isolates were virulent on the Colombian cassava varieties. Such variability in the pathogenicity of Xpm from an international collection without any link between the genetic and pathological patterns was previously reported using different molecular markers [55]. Variations of strain aggressiveness and a lack of correlation between the genetic and pathological patterns were further described at a regional or local scale in different countries revealing dynamic spatiotemporal changes within the Xpm populations [29,30,33,36,38].

Conclusions
The MLVA-12 scheme proposed in this study combined both minisatellite and microsatellite loci that brought complementary information based on their different discriminatory power. Xpm is a monomorphic plant pathogen widely spread in tropical regions. This portable and resolving genotyping tool could be useful to further investigate the global epidemiology of this plant pathogenic bacterium of cassava, whose improved understanding will require additional data.